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Within the framework of a two-state ratchet model, cooperative effects of collective molecular motors on the 
kinetics of adenosine triphosphate (ATP) hydrolysis have been studied in detail. We calculate the dependence 
of the rate of ATP consumption on the applied load. Features different from a single molecular motor have been 
found. We analyze two typical oscillations of collective molecular motors. The relations between the ATP 
concentration and oscillatory characteristics, such as frequency, amplitude, and the duty ratio, have been 


discussed. 


DOI: 10.1103/PhysRevE.69.021912 


Living cells have the ability to move and apply forces. 
Most motions and forces are generated on the molecular 
level by protein molecules that work as motors and are 
driven by a hydrolysis reaction of ATP (adenosine triphos- 
phate) in a situation far from equilibrium [1]. Linear motor 
proteins move along complex periodic and polar structures 
called filaments, obtained by the polymerization of mono- 
mers (actin filament) or dimers (microtubule) [1-6]. Mo- 
lecular motors moving along a one-dimensional polar and 
periodic structure is a nonequilibrium rectification process 
and many different models have been proposed so far to 
describe this process [7-12]. A particularly promising ap- 
proach is due to Vale and Oosawa [13] using the ratchet 
concept introduced by Feynman [14]. Since then, ratchet 
models have been intensively studied. 

Generally, molecular motors are ATPase [1]. The reaction 
kinetics is still an open question. In a recent single molecular 
experiment performed on kinesin (a processive motor) [15], 
the kinetics of ATP hydrolysis has been studied under con- 
trolled external loads. It seems that the kinetics of ATP hy- 
drolysis can still be described by the Michaelis-Menten 
mechanism, which was introduced in 1913 to describe the 
process of enzymatic catalysis [16], although the Michaelis 
constant here is load dependent. The experiment can be well 
explained with a two-state ratchet model. Both the 
Michaelis-like reaction rate of ATP hydrolysis and the veloc- 
ity of molecular motor have been obtained [16]. 

Molecular motors are classified in two groups depending 
on whether they designed to operate in groups(‘‘rowers”’) or 
individually(“porters”) [2]. The collections of motors are 
relevant for experimental situations called “motility assays,” 
for muscle contraction and possibly for vesicle transport in 
cells [2]. Due to collective effects, behaviors very different 
from single individual motor system were found in the col- 
lection of motors, and studied in detail within the framework 
of a two-state ratchet model [2,17,18]. A question arises: 
How does the cooperation affect the kinetics of ATP hydroly- 
sis? In this paper we shall use the two-state ratchet model to 
study the cooperative effects of collective molecular motors 
along a filament on the kinetics of ATP hydrolysis. Features 
of the rate of ATP consumption have been discovered. Oscil- 
latory characteristics, such as the frequency and amplitude, 
as well as the duty ratio have been discussed. 


1063-65 1X/2004/69(2)/021912(5)/$22.50 


69 021912-1 


PACS number(s): 87.10.+e, 05.40.—a 


In the two-state ratchet model [2], a molecular motor is 
represented by a particle. The internal degrees of freedom of 
the motor are simplified to two different states of the particle 
which moves along a one-dimensional x-coordinate axis. 
State 1 is called the “bound” state in which the motor is 
attached to the filament, whereas state 2 is the “free” state in 
which motor is detached from the filament after dissociating 
ATP. The motion of the motor along the x axis at state i(i 
=1,2) is affected by the potential W;(x). For simplicity, 
W,(x) for the “free” state is chosen to be a constant, while 
W(x) for the “bound” state is periodic and with a polarity 
to characterize the filament [see Fig. 1(a)]. The motor transits 
from state 1 to 2 with a rate w,(x), and from state 2 to 1 
with a rate w(x) at position x; both w(x) and w(x) are 
periodic with period / [2,12,16-19]. 

The motor consumes the energy generated by the hydroly- 
sis of ATP to move. The cycle of hydrolytic reaction of ATP 
can be described in a two-step process, as shown in Fig. 2. In 
the first step, the motor gets the energy generated by the 
hydrolytic reaction and tends to transit to state 2 with rate 
a,(x). The transition happens only in the vicinity of the 
minimum of W,(x). In the state labeled “ADPeMeP” the 
motor head is detached from the filament after binding and 


FIG. 1. (a) Ratchet potential W,(x) with U=10kT and W, 
=0.3U=const. (b) @(x) is chosen to be nonzero in the vicinity of 
the minimum of W,(x) with 6//=0.1. (c) Motor molecules are 
connected rigidly to a backbone which is coupled to the environ- 
ment by a spring K; the distance between molecules d is incommen- 
surate with /. 
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FIG. 2. (a) The first step of the ATP hydrolysis take place lo- 
cally. In the vicinity of the minimum of W,, the motor receives 
energy generated by a hydrolytic reaction. (b) The second step is a 
thermal process, i.e, there is not any difference in chemical potential 
for the release of ADP and P. 


dissociating ATP, and the motor may move more or less 
freely along the filament. This state is just the “free” state, or 
state 2. Both state “ATP + M” and state “ADP+M+P” 
refer to state 1. We assume a detailed balance for each step 
[16,20-24]: 


a(x) | Wy(x)— W(x) + Au 1 
a(x) SSE kT , a) 
Pix) = W,(x)— W2(x) 2 
Bo ur a 


where T is the temperature, k is the Boltzmann constant, and 
A u= Mare Lapp  -p is the chemical potential difference. 
The rate of motor transition from state i can be written as 
(x)= a(x) + B(x). 

In order to study collective effects, we model multi-motor 
system as shown in Fig. 1(c). Motor molecules are rigidly 
connected to a common backbone with fixed spacing d 
which is incommensurate with the potential period /. This 
corresponds to the arrangement of motors along their track in 
muscles, and also to the disorder of motors in motility assays 
[2]. Let us denote by P;(&,t) the probability to find a particle 
of state i at position €=x mod / and time t. These densities 
obey P,(&,t)+ P2(€,t) = 1/1 and JIPEN + PE t)dé= 1 
[2,17,18]. The random process of motors along filament sat- 
isfies [2,17,18] 


ðP +tvðP1 = -wP + wP3, (3) 


l 
fen rot | Piadwi-Woae, (4) 


where X is damping coefficient, v=0,X, and f,,, is the ex- 
ternal force acting on the motors. If the system is elastically 
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FIG. 3. For a symmetric system (a/l=0.5) with Xwl?=2kT, 
the velocity (a) and the rate of ATP consumption (b) vs the normal- 
ized concentration of ATP g with different external forces, where 
Ferıl!kT=0.0, —0.02, —0.04, and — 0.08 from top to bottom, re- 
spectively. 


coupled to its environment as shown in Fig. I(c), then f,,, 
= — KX. Note that this model assumes an infinite number of 
motors for simplicity. Recently, computer simulations of fi- 
nite number of motors confirm the validity of this approxi- 
mation [25]. 

Since the hydrolytic reaction of ATP only occurs in the 
first step process, the rate of ATP consumption can be written 
as [16] 


I 
r= [Tarp d-axtdr «le (5) 


It is equal to zero when both Au and f,,, are zero. Because 


TABLE I. The parameters of Michaelis-like relations between the velocity and ATP concentration q at 


different loads. 


Fexıl!kT U max/ ol Ky R Vo/al do 
0.000 0.14887 0.59027 X 104 0.999860 0.11009X 107! 0.5X 10° 
— 0.002 0.12222 0.71157 X 10* 0.999372 0.25394 107! 0.2x 10+ 
— 0.004 0.09287 0.10411 X 10° 0.999560 0.42056 107! 0.5X 104 
— 0.008 0.03893 0.25577 X 10° 0.998476 0.60746 107! 0.3X 10° 
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FIG. 4. For a symmetric system (a//=0.5) with Aoul’=2KT and KI?’=0.02KT. (x;), position vs time; (r;), the rate of ATP consumption 
vs time; (fv;), the external force vs velocity. The subscripts i=1, 2, and 3 correspond to q=10?, 10°, and 10°, respectively. 


the release of ADP and P in the second step [Fig. 2(b)] is just 
a thermal process, we suppose 8>(x) to be a constant (w), 
which is acceptable for W,(x)=const. At the same time, due 
to the “active site” in biology which corresponds the hydro- 
lytic process of ATP, we assume the first step is localized, 
which means a(x) =@6(x), where O(x) is indicated in Fig. 
1(b). Instead of the difference of chemical potential A u, we 
introduce 


_ [ATP] 
7 TADP][P] 


(6) 


which represents the concentrations of ATP normalized with 
respect to their equilibrium value. 

For single molecular motor, both the experiment [15] and 
theoretical analysis of the two-state ratchet model [16] indi- 
cated that the kinetics of ATP hydrolysis can be described by 
the Michaelis-Menten mechanism. The relation between the 
rate of ATP consumption or the velocity of motor and the 
ATP concentrations (q) can simply be given by the Michaelis 
law: 


r= Fm ax4 
Kutg 


(7) 


or 


v = U max] 
Kutg 


tuo, (8) 


where Fmax OF Umax iS the saturation value and Ky is the 
Michaelis constant. For collective molecular motors which 
an fixed external force F,,, is applied, the velocities and the 
rates of ATP consumption versus the normalized concentra- 
tion of ATP q at different loads (f,,,) are shown in Fig. 3. It 
is obvious that both of the velocity and the rate of ATP con- 
sumption cannot be fitted by the Michaelis law. However, at 
a high concentration of ATP, the velocity of motors still sat- 
isfies the modified Michaelis law: 


U masla —4o) 
PUg. 
Kut(g=90) 


v= (9) 
We may call vg and gp the “start velocity ”and the “initial 
concentrations of ATP,” respectively. Table I shows the fit- 
ting of the numerical results with Eqs. (9), and all regression 
coefficients R20.998~1 at different loads. while Ymax de- 
creases when the applied load increases, Ky increases. This 
is similar to the case of a single molecular motor [16]. How- 
ever, the rate of ATP consumption, even at a high concentra- 
tion of ATP, violated the Michaelis law. The rate of ATP 
consumption first increases, then decreases slightly with an 
increasing ATP concentration. This “strange” behavior may 
come from collective effects. The cooperation of molecular 
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motors suppresses the dissipation caused by coupling with 
environment, molecular motors move faster at a lower ATP 
consumption when ATP concentration increases. Figure.3 
also shows that the rate of ATP consumption decreases with 
increasing load. This is in agreement with the single molecu- 
lar motor system [16]. 

Now we consider a system elastically coupled to its envi- 
ronment, as shown in Fig. 1(c). Because of the cooperative 
effects of motors, the system shows spontaneous oscillations 
[2], which explain the oscillatory motions in biological mo- 
tor system, e.g., spontaneous oscillations of single myofibril 
(the contractile unit of muscle cells). We numerically inte- 
grate Eqs. (3), (4), and (5) at fixed damping coefficient X and 
elastic factor K for different concentration of ATP g. The 
establishment of spontaneous oscillations for collective mo- 
lecular motors is shown in Fig. 4 (fv „). From the curve of 
external force (fexr) vs velocity (v) two stages can be seen: 
the transient stage (dash line) and stationary loop stage (solid 
line). Figures 4 show that the motors start at v_ and oscillate 
clockwise along the stationary loop determined by the curve 
of v vs f,xt (dotted line) [2]. This agrees with Jülicher and 
Prost’s prediction [2]. The establishment processes differ at 
low and high ATP concentrations g. At a low ATP concen- 
tration, the oscillation is established gradually, while at a 
high ATP concentration the oscillation is established sud- 
denly. The higher the concentration, the faster the process. 

The rates of ATP consumption (r) satisfy energy conser- 
vation as shown in Fig. 4. In order to maintain the oscilla- 
tion, the system has to offset against the damping force 


( ame — v) by consuming ATP. The energy required 
by the motor per unite time is dE/dt=( —f samping) 


-(dX)/dt=)v?. According to the law of energy conserva- 
tion, the rate of ATP consumption r should be proportional to 
uv’. This is why the frequency of r is twice that of the posi- 
tion. 

Oscillatory characteristics such as frequency, amplitude, 
and the duty ratio of motor at different ATP concentrations q 
are shown in Fig. 5. The amplitude and the time average rate 
of ATP consumption monotonically increase with increasing 
ATP concentration, while the frequency first decreases rap- 
idly and then increases slowly. 

In summary, with the two-state ratchet model, we have 
analyzed the cooperative effects of collective molecular 
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FIG. 5. The frequency (a), amplitude (b), time average of the 
rate of ATP consumption (c), and the duty ratio of the motor (d) vs 
the normalized concentration of ATP q for a symmetric system with 
Nwl?=2kT and KI?=0.02kT. 


motors on the kinetics of ATP hydrolysis. The dependence of 
the rate of ATP consumption on the applied load as a func- 
tion of the ATP concentration has been calculated for cases 
with and without an elastic coupling to the environment. The 
relations between the ATP concentration and oscillatory char- 
acteristics, such as frequency, amplitude, and the duty ratio, 
have been discussed. Features different from those for a 
single molecular motor have been found. 
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